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Abstract 

A set of coupled complex Ginzburg-landau type amplitude equations which op- 
erates near a Hopf- Turing instability boundary is analytically investigated to show 
localized oscillatory patterns. The spatial structure of those patterns are the same 
as quantum mechanical harmonic oscillator stationary states and can have even or 
odd symmetry depending on the order of the state. It has been seen that the under- 
lying Turing state plays a major role in the selection of the order of such solutions. 
PACS number (s): 87.10.+e, 47.70.Fw 

Localized spatio-temporal structures play a very important role in the long time dy- 
namics of many chemical and hydrodynamic systems. Investigation of local instabilities 
is an important field [1-8], specially the so called topological defects in chemical and bi- 
ological systems [9-15] which are believed to play a role in the transition from phase to 
defect chaos regime. Interaction of local and global modes give rise to very complicated 
large scale dynamics [16-19]. It is still an unsolved problem as to how dissipation and 
nonlinear feedback balance each other to produce localized oscillations (oscillons) and in 
that context some field theoretic studies [20-24] are worth mentioning. So, it makes a 
very rich and broad field - the interaction of very small and large scales- in the context 
of dynamical systems. In an extended system, the formation of spatio-temporal defects, 
which are generally localized in nature on a stationary or uniformly moving frame, are 
generally believed to be generated as a result of local and global symmetry mismatch 
which can arise from boundary conditions, multi-stability of nonlinear states etc. On the 
contrary, intrinsic local instabilities of a system with its own symmetry can act as a sym- 
metry breaking agent of the global state its embedded in and are very important in this 
respect. By intrinsic we mean those which have not been generated by multi-stability of 
global states such as a droplet of one state pinned by suitable front structures within an 
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extended one. Rather, these are local instabilities to the global states. Such instabilities 
can grow up locally and be connected to a global state by fronts or pulses whose symmetry 
reflects that of the localized core. Thus, its about thinking the other way where basic 
symmetry of the localized core that determines the symmetry of the large scale structures. 

In the present work we are going to show a class of new local instabilities which are 
like harmonic oscillator states and can arise in oscillatory as well as steady form. The 
most interesting thing about this class of local structures is their spatial symmetry which 
can be even and odd depending on the order of the state. Since these local structures 
are Hermite polynomial in their spatial profile, they can be symmetric or antisymmetric 
depending upon the order of the polynomial. In the present work it has been shown that 
instabilities of order greater than zero can grow depending upon the situation prevails. 
Such symmetric and antisymmetric local structures are also shown to be a steady state 
of the system with relatively larger spatial extent. We will show that for such a local 
instability to grow, a coupling of the Turing and the Hopf modes are essential when 
the bifurcation is supercritical. A pure Hopf environment can also produce oscillatory 
as well as steady (more restricted) form of such localized structures. Although, in the 
following we mainly work in a region of phase space close to a Hopf- Turing instability 
boundary but the result can be generally applied to any situation governed by a complex 
Ginzburg-Landau equation having the real part of the coefficient of the diffusive term to 
be negative. In a purely Turing region, it will be shown that, the instability will appear 
in a non-oscillatory form and a coupling to Hopf mode is essential for it to oscillate. 

The slowly varying amplitudes of Turing and Hopf instabilities being denoted by T 
and H the coupled amplitude equations in one dimensions are given by 

ea t T + b t -^ - c t \T\ 2 T - d t \H\ 2 T 

ea h H + b h ^r-c h \H\ 2 H-d h \T\ 2 H (1) 



dT 

~dt 
dH 

~dt 



Where all the coefficients with a suffix h are generally complex and those with suffix t are 
real. This system has a steady state of the form T = T e lkx and H = H e tult where the 
wave number and the oscillation frequency are given by 

2 _ ea t — Q | To | 2 — d t \H \ 2 

and, 

^ = £ a hi — Chi\H \ — d hi \T \ 

where the suffix hi represents the imaginary part of the complex coefficients. Let us see 
what happens to local perturbations which has a Gaussian profile to these steady states. 
To see that, we put the perturbations of the form STe~ x ^ 2b and 5He~ x ^ 2b on the Turing 
and Hopf states respectively to get 

- d t H T Cos(ujt)e ikx 5H 



85T 



= KrfT + bt 



d 2 5T 2x d5T 



dx 2 



b dx 



2 



dSH 
~df 



K 2 5H + b h 



d 2 5H 2x dSH 



dx 



b dx 



- d h H T Cos{kx)e iwt ST 



(2) 



where 



K x = ea t -^ + b t ^-c t T 2 -d t H 2 
b b 2 

K 2 = ea h - y + b h ^ - c h H 2 - d h T 2 



(3) 
(4) 



The region of investigation is such that x 2 /b 2 is very small and we will neglect it. Let us 
re-scale the space co-ordinate such that x — and thus we are looking at a somewhat 
bigger scale with b > 1 and then make some reasonable approximations. The second 
approximation is - region of local structures are so small compared to the long wavelength 
of the Turing order k we will replace Cos(kx) and e lkx by unity. This approximation is 
justified because comparatively large scale modes are supported by amplitude equations. 
With these simplifying assumptions and rescaling the linearized model looks like 



d5T 
05H 



b 

= K 2 5H+% 
b 



d 2 5T _ 
dx 2 
d 2 5H 



2x- 



d5T 



2x 



dx 
dSH 



dx 



d t H T Cos{ut)SH 
- d h H T e wt 5T 



(5) 



where now K 1 and K 2 are devoid of terms containing x 2 /b 2 . To separate space and time 
let us take the solution in the form 



5T 
5H 



T(t)H n (x) 
H(t)H n (x) 



(6) 



where H n (x) is the Hermite polynomial of nth order. Putting this in the above linearized 
Eq.5 and dividing them by ST and 5H respectively we get 



1 dT 
Tifj'dt 
1 dH 
H{t) dt 



= K,+ 
= K 2 + 



bH n (x) 
bH n (x) 



d 2 H n (x) _ 2x dH n (x) 



dx z 
d 2 H n (x) 
dx 2 



2x 



dx 
dH n (x) 
dx 



dh,H T 
H(t) 



Cos{ut)H{t) 
^Tit) 



(7) 



Each of the above two equations now can be separated in the space and time parts. After 
we separate the space and time part of above equations, we see that the space parts are 
clearly decoupled and the coupling exists only in the time part, the decoupling of the 
space part is very important in view of the fact that we can now set them equal to some 



3 



arbitrary constants and let them be — 2ny and — 2n^ respectively for the Turing and 
Hopf case where n is an integer. The space parts would now look like 



bH n (x) 

h 
bHJx) 



d 2 H n (x) 

dx 2 
d 2 H n (x) 

dx 2 



2x 
2x 



dHnjx) 

dx 
dH n (x) 

dx 



b 

-%2n 
b 



(8) 



Thus, Eq.8 makes it visibly clear that in each part the equation for Hermite polynomial 
solution is satisfied and there is a spatial instability in the form of a harmonic oscillator 
state e~ x l 2 H n (x) (for 6=1) and the condition for its growth is to be determined. 
The temporal parts give coupled linear equation as 



dT(t) 

~dJ~ 

dH(t) 




T(t) - d t T H Cos(ut)H(t) 
H{t) - d h T H T(t)e tuJt 



(9) 



Now we would like to see whether the spatial solutions of the local structures in the form 
of various bound states of a quantum mechanical harmonic oscillator grows or not. We 
know that K\ = b t k 2 — y and K 2 = —^- + iuj where the second relation comes from 
the steady state condition of the Hopf mode. With K x and K 2 replaced by appropriate 
expressions Eq.9 looks like 




b t (k 2 -^) -d t T H Cos(ujt) 
-d h T Q H Q e™ -b h ^ + icu 



(T{t) 
\H{t) 



(10) 



Thus, with the coupling absent a steady instability will grow in the amplitude of the 
Turing mode under the condition 



k 2 > 



2n+ 1 



(11) 



where the wave number of the preexisting rolls decides on the upper limit of the order of 
localized structure given the fixed spread of it. Although we see here that an instability 
of lowest order n = has maximum linear growth, but that may not be the only selection 
criterion for the order of the instability. Other factors and nonlinearity might play in 
selecting a nonzero order when they come in. When b t is negative the situation reverses 
and the wave number k gives here the lower bound to the order n and we can safely say 
that in the case of growth of an instability in the Turing part when 



k < 



2n + 1 



(12) 



In such a situation the nonlinearity must play a role in putting the upper bound to the 
order of the instability since linear growth is proportional to n. Since the bifurcation 
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is sub critical with negative b t which is more effective in producing nonzero order of the 
instability, the transition can show hysteresis. Now, what the coupling of Turing and Hopf 
mode does is making the steady instability developed in pure Turing situation oscillate 
and obviously producing some sort of frequency entrainment with the underlying Hopf 
mode following the standard theory of forced linear systems. 

In the same decoupled situation with the Hopf part of the equations we see that it also 
admits growth of such an oscillatory instability separately but only when the real part of 
bh is negative. With the real part of bh equal to zero the growth of a steady instability is 
marginal with an order n satisfying the relation y- — 2n + 1. This is a rather stringent 
condition which allows a discreet set of u of the preexisting Hopf mode to produce steady 
local structures. Thus, a continuous band of uo is effective in making grow an oscillatory 
instability with an oscillation frequency 

cu=(2n+l)^-cu (13) 
b 

and it would be more commonly observed in situations governed by a single one di- 
mensional complex Ginzburg-Landau equation. No coupling with the Turing mode is 
necessary here to produce such an oscillatory instability to grow but on the other hand 
the coupling to hopf mode is essential for an instability to oscillate which originates from 
the underlying Turing state. Also in the Hopf part the instability of higher order has a 
higher growth rate and one would expect to see some with nonzero n and the other factors 
have to play a role on the selection of the upper bound to the order. 

The important fact about these instabilities is that they are intrinsic in nature and can 
appear separately in Turing or Hopf environment. At the same time, a mixed Hopf- Turing 
phase is also locally unstable to such instabilities. In a large b limit such steady harmonic 
oscillator state like local structures are a solution of the system when one could neglect 
the nonlinear terms away from the origin where e~ x l h is very small and thus, a region 
away from the origin but not too far will be appropriate for the appearance of these steady 
structures. In such a situation the order of the steady state is given by ea t = (2n+l)y and 
similar relation will come from the Hopf part by separately equating real and imaginary 
coefficients. Thus, the Hopf part is capable of producing bi stability of such steady local 
structures when the real and imaginary parts differ in magnitude. However, a detailed 
treatment of such symmetric and antisymmetric local instabilities giving rise to moving 
pulse or front structures or surviving in the presence of inhomogeneous oscillatory global 
environment can prove to be very interesting due to the clash of symmetries. Its the 
spatial symmetry of these structures which can be odd and would definitely make the 
coexistence of them in a spatially symmetric global state an interesting thing to study. 
When oscillates, such antisymmetric localized structures can also influence the dynamics 
of its surroundings in a nontrivial way and further work in the field of interaction of these 
localized structures with its environment can reveal such effects. 
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